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Abstract. An exact solution is derived for the wave function of an electron in a 
i-Ch \ semiconductor quantum wire with spin-orbit interaction and driven by external time 

dependent harmonic confining potential. The formalism allows analytical expressions 
for various quantities to be derived, such as spin and pseudo-spin rotations, energy 
and occupation probabilities for excited states. It is demonstrated how perfect spin 
and pseudo-spin flips can be achieved at high frequencies of order u), the confining 
potential level spacing. By an appropriately chosen driving term, spin manipulation 
can be exactly performed far into the non-adiabatic regime. Implications for spin- 
polarised emission and spin-dependent transport are also discussed. 
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One of the first proposals for the solid state realization of quantum information 
processing and computation is based on the use of the spin degree of freedom of electrons, 
confined in quantum dots (QD), which realize the qubit pQ, and several universal 
quantum computation propositions [21 El H] have been deduced from this basic idea. 

Experimental progress in measurements and coherent manipulation of electrons in 
QD systems has been immense in the past years. El El IHl El EDI dH (121 ESI H31 HSl dSl [12] 
Coherent single electron spin rotations have been realized using oscillating magnetic [5] 
and oscillating electric fields. j6j [7] The limitation of these schemes is due to long single 
spin rotation times, of the order of tens of ns. By exploiting ultrafast laser control 
and optical pumping, single spin control was demonstrated on a ps timescale.jEl E] 
All-electrical two-qubit gates composed of single-spin rotations and interdot exchange 
in a double QD were also demonstrated. [10] By using surface acoustic waves, single 
electrons have been transfered between quantum dots separated by a few /im [HI [12], 
where McNeal et. al. have reliably relocated a single electron back and forth between the 
dots multiple times up to a cumulative distance of 0.25 mm.[T2"] Controlled production 
of entangled qubit pairs in quantum dots and surface acoustic waves was predicted 
theoretically. [181 [T9| [20] Coherent time evolution of single electron charge dynamics on a 
time scale of a few picoseconds was measured. [13] Tunable non-adiabatic single electron 
excitations were observed recently, [H] where rapid modulation of confining potential 
causes transitions to excited orbital states. 

Due to the lack of inversion symmetry in semiconducting materials, the electron 
spin is coupled to the orbital motion. [2TJ [22] Rather than treating this interaction as 
a source of decoherence, it can be harnessed to coherently manipulate spin-orbital 
qubits.jl EH EU EH EE E3 EE] The rotation of a spin-orbital qubit in a moving 
quantum dot has been studied analytically in the adiabatic limit. Single qubit 
manipulation, spin-flip for example, can easily be performed by adiabatic spatial QD 
translation for a distance of the order of the spin-orbit length. [2H EHJ EHJ E7J However, 
the condition of adiabaticity sets a severe upper limit for the speed of single qubit 
operations. Although the analysis by means of a full numerical simulation is possible, 
it gives results only for particular cases and the optimization and elimination of errors 
is not straightforward. [27] 

In this paper we extend analytic results for spin manipulation by adiabatic 
translation of a QD [211 [251 EE] to include the non-adiabatic regime, thereby opening the 
possibility of much faster spin-qubit rotations. In particular, we analytically solve the 
one-dimensional (ID), non-adiabatically driven QD with harmonic confining potential 
in the presence of spin-orbit interaction. Although it is demonstrated how a general 
solution can be constructed, we also present explicitly four typical driving schemes with 
the focus on the non-adiabatic regime. The system considered here is similar to the one 
presented by Flindt et. al.,[23\ where the authors performed a spin-flip by displacement 
of a QD in a wire. Such a system has even been realized experimentally, EH [291 ETJ1 EI] 
with the ability to control a single electron, in InAs [30] and Ge [31] quantum wires. 
More recently, spin-orbital qubits in double QD have also been demonstrated in InSb 
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Figure 1. (a): Schematic view of a semiconductor quantum wire deposited on 
gate electrodes which provide a time dependent harmonic potential. Such systems 
have been realized experimentally [7J [28] [29]. (b): Position of the harmonic 
potential minimum. After driving time T the QD stops at displacement £(T). Black 
line corresponds to the fastest unidirectional QD translation, with minimum spin-flip 
time To/2, where To = 2-k/lo and u> is the confining potential level spacing. Besides 
instantaneous displacement, the other drivings considered are: linear ramp with = 
£(T)t/T (dashed blue) and sinusoidal with = £(T) [t/T± sin(2 nt/T)/ (2 tt)] (dot- 
dashed green and solid red), giving spin- flip times of T , T /\^2 and 2 To, respectively. 



quantum wires, j 

The ID Hamiltonian we consider is given by 

H(t) = ^ + ^f- [x - +P-A»a y -(3 a x ), (1) 

where an electron with an effective mass m* is confined in a harmonic trap with frequency 
oj, centered around the time dependent position of harmonic potential minimum £(£), a 
and (3 are the Rashba (21] and Dresselhaus [22] spin-orbit couplings, respectively, and 
we use units with h~ = 1 throughout the paper. 
First we apply the unitary transformation 

U\t) = e- s V\t), (2) 

which consists of a canonical transformation that eliminates the spin-orbital part of the 
Hamiltonian, [32] 



S = i m* x (aa y — (3 a x ), (3) 
and of the part which eliminates the time dependent part of the Hamiltonian, [33] 
V\t) = e 1 m * ±c ^ x e ~ iXc ^P x ^ 

m = - f mat, (4) 
L( t ) =^x c (tf-^f- [ Xc (tf-m 2 ]- 

Here x c (t) is the solution of the equation of motion for a classical particle in a driven 
harmonic oscillator potential, 

x c + u 2 x c = w 2 £(t), (5) 
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for x c (0) = and x c (0) = customarily given by x c (t) = u f Q sm[co(t — t')] £(t') dt'. 

The time dependent Schrodinger equation after the unitary transformation 
U (t) H (tyU^ (t) — > H takes the form of a static harmonic oscillator 

iM= I^ + ^J + eJ*, (6) 
vim* 1 J 

with an energy shift E so = —m* (a 2 + (3 2 )/2. The exact eigenstates of the original 

Hamiltonian ([TJ), ty ns (x,t), are obtained directly from the eigenstates of if) n (x,t), 

via the unitary transformation ([2]), z.e., 

V ns (x,t) = e-^+^U^^ix) Xs = 

= e i ^^i/j n [x-x c (t)}e- S Xsi (7) 

where <p n (x, t) = —{E so + u n )t + $(£) — m*x c (t) x is a phase factor, w n = (n + l/2)u> is 
the eigenenergy of the n-th eigenfunction of the static harmonic oscillator ip n , and Xs a 
spinor with spin s. 

From the exact time-dependent solution, (JJJ), we see that the effect of non-adiabatic 
motion is to displace the harmonic oscillator ground-state gaussian wavepacket away 
from the quantum-dot minimum at £(t) to a new position x c (t), given by the solution 
of the classical oscillator, (JHJ). It is clear that this will excite the electron in the QD and 
explicit calculation gives 

S n (t) = {y ns \H(t)\y ns ) = 

_ m* x c (t) 2 m* u 2 , . . 1 _. ,.o , , 

= E so + u n + + [x c (t) - £(t)} \\ (8) 

which is in fact the classical expression for the harmonic oscillator energy, apart from 
a constant shift of E so + u n . To simplify analysis we consider quantum wires grown 
in the [111] direction, for which /3 = 0. For this case the mean spin rotates around 
the y-axis by an angle 2 x c (t)/ \ so , where \ so = (m* a)" 1 , following immediately from 
equation ((7j), from which the expectation values of the spin components may also be 
evaluated explicitly. For example, if the initial state is the pseudo-spin-up Kramers 
state, ty(tf(x, 0) = e ~ iX / XsorTy ipo(x) xt> the spin expectation values are 

S,(t) = \{a z ) = l - C os[2x c {t)/X so }e- 2a2 / x "% 

S x {t) = l -(a x ) = l -sm[2x c {t)/\ so ]e- 2 ° 2 l x2 s % (9) 

where the oscillatory factors are consistent with the mean spin rotation and the decaying 
exponential factors account for the spread in the wavefunction, which is given by 
a = (m* w)~ 1 / 2 . Due to the precession of spin around the y-axis S y is constant and 
zero for our choice of the initial state. 

In figure DJb) we show specific choices of driving term, for cases when the 
transit time, T, is chosen to give a perfect spin-flip, for which the dot displacement is 
ir\ so /2. Detailed results for these choices and general T are given below. Here we note 
that in general there will be residual oscillations of the electron distribution after the dot 
stops, leading to uncertainty of order 2a/X so in the rotated spin rotation angle , where 
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a is the amplitude of the residual orbital oscillations. Such residual oscillations of the 
electron wave packet after the dot has stopped will also give rise to orbital decoherence 
through phonon emission and subsequent spin decoherence via the spin-orbit interaction. 
However, it is clear that £(£) can always be chosen such that these residual oscillations 
are totally suppressed with the wavefunction returning to the centre of the QD and 
remaining there, i.e., x c (t) = £(T) and x c (t) = for all t > T. This is a very physical 
result which can be understood in purely classical terms from equation (JSJ). The final 
energy is given by £ fina i = E so + u n + m* u 2 a 2 /2. 

A further interesting general observation is that under forced oscillations in which 
the quantum dot is translated back and forth non-adiabatically in one dimension, energy 
will be transferred to the electron wave packet in the dot and under resonant conditions 
the amplitude of the displacement of the classical oscillator (and hence the position 
of the wavefunction relative to the potential minimum) will increase, as will the mean 
energy of the electron given by equation (151). Accompanying these oscillations will be 
a corresponding rotation of the spin which will also have oscillatory components, due 
to the oscillation of the dot and oscillation of the wavefunction relative to the dot. It 
would be interesting to observe effects of these oscillations experimentally, e.g. in the 
polarisation distribution of radiative emission when the excited electron in the quantum 
dot relaxes or by using weak repeated measurement of current. [T3l [T3] 

In materials with spin-orbit interaction, besides the real spin, one can study also 
the pseudo-spin, defined by 

^ = 2 E(l^«t)(^t|-I*n;)(*n;|), 

n=0 

t = - ]Cd*«tX^I + I^X^ntl), (10) 
n=0 

where ^ ns (x, t) are the eigenstates of the potential at a particular instant t, 

V ns (x, t) = e~ im * a ^ i> n \x - £(*)] Xs- (11) 

For the initial Kramers pseudo-spin-up state, a straightforward calculation leads to a 
simple exact expression for the expectation value, 

W) = (%) = ~cos[2£(f)/A so ], 

T x (t) = (%} = ^sin[2 ^t)/ X so ], (12) 

and T y = 0. Pseudo-spin is dependent only on the position of the QD, therefore slow 
adiabatic or arbitrary fast non-adiabatic motion lead to identical pseudo-spin values. 
Furthermore, the pseudo-spin is independent of the spread in orbital wavefunction, 
being a composite of spin and orbital degrees of freedom. In the derivation of this result 
we express the coefficients in the instantaneous basis given by fill I) . 

^(O^^l^^e-^+^^lxle-^WA-^xt^), (13) 
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/„(*) = J dxe- im ' ± ^ x r n [x-m}M^-^(t)} = 

= e -x c {t) 2 /(2a 2 w 2 )+im* x c {t)x c {t) ^ n |Q\ d) 

where (n\0)d = f duij}*(u) ipoiu — d) is the overlap of the n-th eigenfunction of 
the quantum harmonic oscillator and the ground state displaced by the distance 
d = Xc (t) - £(t) -ix c (t)/u, and (n\0) d = [d/(V2a)] n /VnJ. exp[-rf 2 /(4a 2 )]. Note that 
this result follows due to the independence of the spin and the orbital part. Thus the 
probability of occupation of the n-th manifold is exactly determined solely by classical 
coordinates, P n (t) = J2 S \c ns (t)\ 2 = \I n (t)\ 2 . 

From equation (Tl3l) we see that the exact wavefunction at time t may be written, 

ns 

= e -* [(*«-*"/»>*-*('>] [cos(£(t)/A so ) tf t + sin(£(t)/A S0 ) (14) 

where = J2 n In{t) *&ns(x, t). In this form we see explicitly the rotation of the pseudo- 
spin on the Bloch sphere and also the time-dependence of the orbital part of pseudo-spin 
states \& s . 

Pseudo-spin qubits in systems with strong spin-orbit interaction are usually 
restricted to the manifold spanned by the ground state Kramers doublet, as studied in 
references [7J [231 EH EHJ EU [27] . Pseudo-spin fTTOT) is in this case limited only to the n = 
manifold. For the initial Kramers pseudo-spin- up ground state ^/ot( x 5 0), the expectation 

l 2 i 2 

value of pseudo-spin is related to the probabilities Pot = l c otl an d Poi — \ c oi\ that the 
system is in the ground state with pseudo-spin up or down, respectively. 

Direct evaluation gives T z °(t) = |(P f ~ Pol), where P f = c os 2 (^/A so ) P , P oi = 
sin 2 (£/A so ) P , which leads to 

7?{t)=Ti(t)P (t), i = x,y,z, (15) 

and the probability for the ground state manifold occupation is given by Po(t) = 
exp[— (x c — ^) 2 /(2<t 2 ) — x 2 c / {2 a 2 uj 2 )}. Note that in the adiabatic limit the moving 
electron remains at the potential minimum, i.e., x c (t) — > £(£) and Po — > 1 and the 
pseudo-spin coincides - up to an exponential factor due to a wavefunction spread - with 
the real spin, given by Qj. 

To be specific, we examine four examples of the driving term £(£): sudden 
displacement in two steps, a constant velocity of the QD displacement, and two different 
types of smooth driving with sinusoidal (AC) type of velocity, as shown in figure d^b). In 
all cases the QD stops moving after the transit time T. Then, the position expectation 
value oscillates with the amplitudes of residual oscillations a, as follows from the classical 
solution x c (t) of the equation of the motion (jHJ). For example, 

Aso Tq 



2 T 



T J 
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for constant velocity (a c ) case and smooth sinusoidal case (a s ), as introduced in figure 
CD (solid red). 

The shortest time with no residual oscillations for unidirectional movement of the 
dot, T min , is achieved by instantaneous displacement of the dot for 7rA so /4, followed 
by a waiting time To/2, followed by another instantaneous displacement for 7rA so /4. 
For any other unidirectional driving the minimum time T min rises. In the limit 
T —7- all unidirectional movements of the dot lead to oscillations with the amplitude 
£(T) = it \ so /2. In figure [2(a) residual oscillations are presented as a function of driving 
time T. In all cases the residual oscillations vanish when some resonant condition for 
T/Tq is fulfilled. The oscillations are proportional to the spin-orbit length X so and have 
a simple dependence on T an d T. As expected, smooth sinusoidal driving exhibits 
significantly less residual oscillations compared to other forced oscillations considered. 
As demonstrated here, by using the solution to a classical driven harmonic oscillator we 
gain insight into a driven QD in a semiconducting nanowire. 
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Figure 2. (a): 



Amplitude of residual oscillations after the QD stops, a = 
2 divided by the QD displacement £(T) = ir\ so /2 for 



drivings as in figure 1(b): Fastest unidirectional driving (full black), constant velocity 
(full red) and the two sinusoidal drivings, a broken (dot-dashed green) and a smooth 
one (dashed blue), (b): Probability for the ground state doublet occupation after 
the QD stopes. The material parameters that we use are given by m* = 0.015 m e 
and Xso = 150 nm, where m e is electron mass, which are typical for InSb, [28], the 
confining potential level spacing w = 10 meV and a = 0.1 £(T). 



In figure [3]^a) the ^-component of spin after one cycle is presented for linear ramp 
and smooth sinusoidal driving for T/T = 2 and 5. For this spin-flip case the quantum 
dot is displaced a distance £(T) = ir A so /2. Note that the magnitude of S z does not 
quite reach | at the beginning and end of the spin-flip due to the spread in the gaussian 
wavefunction, which gives S z (0) = \ exp (— 2 cr 2 /\ 2 so ) from equation Qj. Note also the 
weak oscillation superimposed on the cosine curve which, again, can be understood 
directly from the exact solution, equation OH]), which shows a perfect cosine variation 
when plotted against x c rather than £. For x c — £ small we may expand to lowest order 
giving cos(2x c /A so ) = cos(2£/A so ) — 2 (x c — £) sin(2£/A so )/A so which is the form shown 
in the figure, since x c — £ has n = T/T oscillations. These oscillations are both positive 




Figure 3. (a) z-componcnt of spin expectation value as a function of QD displacement 
in non-adiabatic regime for two types of driving. Full black line represents the adiabatic 
limit (T — > oo). Note the oscillations superimposed on the cosine curve giving rise to 
both positive and negative fluctuations in mean spin relative to the adiabatic limit, 
(b) Ground-doublet pseudo-spin expectation, 7j\ for parameters as in (a). Note that 
T 2 °(0) = -T 2 °(T) = \ precisely and that the effect of non-adiabaticity leads to greater 
deviation from cosine (compared with real spin) at intermediate times, which is always 
suppressed due to breakdown of the ground-doublet approximation, (c): Typical 
probability for occupation of the ground state doublet during the driving (sinusoidal), 
Po(t), in the non-adiabatic regime. The arrow indicates the direction of increasing T 
(T/Tq=2, 2.5, 3.5, 4.5, 6.5); bullets indicate minimum values Pomin- (d): -Fbmin &s a 
function of T/Tq; bullets as in (c). All parameters as in figure EJb). 
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and negative with respect to the adiabatic limit (full black line) which shows a perfect 
cosine (also with amplitude attenuated by the gaussian factor in equation (jUJ)) but with 
no superimposed oscillation since = x c (t) for all t. 

In figure [3](b) the ground- doublet pseudo-spin expectation 7^° is presented for the 
same conditions as in (a). The adiabatic limit (full black line) is now a perfect cosine 
with amplitude | and this is also the result for full pseudo-spin T z for any T consistent 
with spin- flip (see equations (fTUl) and (I12p ). On the other hand, we see that there are 
significant deviations for non-adiabatic motion in the ground-doublet approximation, 
signaling its breakdown due to suppressed probability of ground-doublet occupation 
and corresponding suppression of 7^° relative to the adiabatic limit. These deviations 
do not occur, however, at the beginning and end of the spin-flip for which the function 

have been chosen so that the electron dot starts and ends in its ground- doublet 
located in the centre of the dot. Nevertheless, we see from the figure ©a) and[3^b)) 
that any error in the total displacement will give spin-flip errors which are significantly 
smaller for the smooth sinusoidal choice of as expected. 

For all the different driving functions the optimum spin (or pseudo-spin) 

transformations, which leave the electron in its lowest energy state after transit of the 
dot, may always be determined. Such states depend only on the distance the dot travels, 
independent of the degree of non-adiabaticity, provided the motion is coherent. However, 
the states of the electron in the dot at intermediate times/distances do depend strongly 
on the degree of non-adiabaticity and a convenient measure of this is the deviation from 
unity of the probability for the ground state, Pq. For t > T this probability is simply 
determined by the residual oscillations, Po(T) = e~ a2 ^ 2a ' 2 \ as shown in figure [2(b) 
for driving functions as in figure EJ^a). In figure [Ujc) the time dependence of Po(t) for 
different values of T/T is presented for the smooth sinusoidal case: for large T, P ~ 1, 
while for short times the minimum of P = Pomm around t ~ T/2 signals the increase in 
probability for exciting higher energy states. The transition from highly non-adiabatic 
to the adiabatic regime is clearly seen also from the behaviour of the minimum value of 
Po, plotted in figure EJ^d), which at high T monotonically reaches the unitary limit. 

To conclude, we have presented exact solutions of an electron in a driven QD in a ID 
quantum wire with spin-orbit interaction. The solutions can be expressed analytically 
for a broad class of driving schemes. The essence of the full quantum mechanical solution 
is in the solution for the corresponding classical harmonic pendulum and by applying 
the appropriate canonical transformation. Exact analytical formulae for spin, pseudo- 
spin, average energy and probabilities for the occupation of excited states are given 
in terms of the coordinates of the QD and the position of the corresponding classical 
pendulum. To be specific, we mainly focused on two typical driving regimes: a constant 
velocity and a smooth sinusoidal driving - in both cases in the full range from adiabatic 
to non-adiabatic driving. We demonstrated that perfect spin-flip can be realized even 
in the far non-adiabatic regime, providing driving which by the virtue of the formalism 
can easily be appropriately tuned. The solution given here may be used to guide more 
complicated simulations for the design of realistic gated device structures. 
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Two extreme non-adiabatic regimes were highlighted. The first is when the forced 
oscillator is chosen such that after the transit time of the quantum dot the wavefunction 
of the electron is at the bottom of the QD potential well and remains there. The main 
advantage of this choice is that the associated spin rotation (single-qubit transformation) 
may be determined with maximum precision and does not subsequently change with 
time. The other extreme regime is exactly the opposite in which the forced oscillations 
have maximum amplitude after the transit time (resonance). Indeed, we can increase 
this amplitude indefinitely using a periodic driving term for the dot motion when it 
is resonant with the harmonic oscillations of the electron in the QD. The amplitude 
of these oscillations will, of course, be ultimately limited by dissipative processes, in 
particular phonon and radiative emission, though in clean systems coherent oscillations 
should survive to large amplitudes. What is special about such systems is the intimate 
connection between motion of the dot, the electron probability distribution relative to 
the dot and spin rotation. It would be interesting challenge to observe the effects of 
such charge-spin oscillations experimentally, for example in spin-polarised emission or 
spin-dependent transport. 

Acknowledgments 

The authors thank T Rejec for discussions and acknowledge support from the Slovenian 
Research Agency under contract Pl-0044 and the EU grant NanoCTM. 

References 

[1] Loss D and DiVincenzo D P 1998 Phys. Rev. A 57 120 
[2] Levy J 2002 Phys. Rev. Lett. 89 147902 

[3] DiVincenzo D P, Bacon D, Kempe J, Burkard G and Whaley K B 2000 Nature 408 339 
[4] Stepanenko D and Bonesteel N E 2004 Phys. Rev. Lett. 93 140501 

[5] Koppens F H L, Buizert C, Tielrooij K J, Vink I T, Nowack K C, Meunier T, Kouwenhoven L P 

and Vandersypen LMK 2006 Nature 44 2 766 
[6] Nowack K C, Koppens F H L, Nazarov Y V and Vandersypen LMK 2007 Science 318 1430 
[7] Nadj-Perge S, Frolov S M, Bakkcrs E P A M and Kouwenhoven L P 2010 Nature 468 1084 
[8] Berezovsky J, Mikkelscn M H, Stoltz N G, Coldren L A and Awschalom D D 2008 Science 320 

349 

[9] Press D, Ladd T D, Zhang B and Yamamoto Y 2008 Nature 456 218 
[10] Brunner R, Shin Y S, Obata T, Pioro-Ladricrc M, Kubo T, Yoshida K, Taniyama T, Tokura Y 

and Tarucha S 2011 Phys. Rev. Lett. 107 146801 
[11] Hcrmclin S, Takada S, Yamamoto M, Tarucha S, Wicck A D, Saminadayar L, Baucrle C and 

Meunier T 2011 Nature 477 435 
[12] McNeil R P G, Kataoka M, Ford C J B, Barnes C H W, Anderson D, Jones G A C, Farrcr I and 

Ritchie D A 2011 Nature 477 439 
[13] Kataoka M, Astlcy M R, Thorn A L, Oi D K L, Barnes C H W, Ford C J B, Anderson D, Jones 

G A C, Farrcr I, Ritchie D A and Pepper M 2009 Phys. Rev. Lett. 102 156801 
[14] Kataoka M, Fletcher J D, See P, Giblin S P, Janssen T J B M, Griffiths J P, Jones G A C, Farrcr 

I and Ritchie D A 2011 Phys. Rev. Lett. 106 126801 



Non-adiabatically driven electron in quantum wire with spin-orbit interaction 



11 



[15] Petta J R, Johnson A C, Taylor J M, Laird E A, Yacoby A, Lukin M D, Marcus C M, Hanson 

M P and Gossard A C 2005 Science 309 2180 
[16] Hanson R, Kouwenhovcn L P, Petta J R, Tarucha S and Vandersypcn L M K 2007 Rev. Mod. 

Phys. 79 1217 

[17] Bennett A J, Poolcy M A, Stevenson R M, Ward M B, Patcl R B, dc la Giroday A B, Skold N, 

Farrer I, Nicoll C A, Ritchie D A and Shields A J 2010 Nature Phys. 947 
[18] Giavaras G, Jefferson J H, Ramsak A, Spillcr T P and Lambert C J 2006 Phys. Rev. B 74 195341 
[19] Jefferson J H, Ramsak A and Rejcc T 2006 Europhys. Lett. 74 764 
[20] Rejcc T, Ramsak A and Jefferson J H 2000 J. Phys.: Condens. Matter 12 L233 
[21] Bychkov Y A and Rashba E I 1984 J. Phys. C 17 6039 
[22] Dresselhaus G 1955 Phys. Rev. 100(2) 580 

[23] Flindt C, S0rcnscn A S and Flcnsbcrg K 2006 Phys. Rev. Lett. 97 240501 

[24] Coish W A, Golovach V N, Egucs J C and Loss D 2006 Phys. Stat. Sol. (b) 243 3658 

[25] San- Jose P, Scharfenberger B, Schon G, Shnirman A and Zarand G 2008 Phys. Rev. B 77 045305 

[26] Golovach V N, Borhani M and Loss D 2010 Phys. Rev. A 81 022315 

[27] Bednarek S and Szafran B 2008 Phys. Rev. Lett. 101 216805 

[28] Nadj-Perge S, Pribiag V S, van den Berg J W G, Zuo K, Plissard S R, Bakkers E P A M, Frolov 

S M and Kouwenhoven L P 2012 Phys Rev Lett 108 166801 
[29] Fasth C, Fuhrer A, Bjdrk M T and Samuelson L 2005 Nano Lett. 5 1487 
[30] Fasth C, Fuhrer A, Samuelson L, Golovach V N and Loss D 2007 Phys. Rev. Lett. 98 266801 
[31] Shin S K, Huang S, Fukata N and Ishibashi K 2012 Appl. Phys. Lett. 100 073103 
[32] Khaetskii A V and Nazarov Y V 2000 Phys. Rev. B 61 12639 

[33] teer Haar D 1964 Selected problems in quantum mechanics (Infosearch Limited, London) 



